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The Fermi satellite has been reporting the detailed temporal properties of gamma-ray bursts (GRBs) in an 
extremely broad spectral range, 8 keV - 300 GeV, in particular, the unexpected delays of the GeV emission 
onsets behind the MeV emission of some GRBs. We focus on GRB 080916C, one of the Fermi-LAT GRBs for 
which the data of the delayed high-energy emission are quite extensive, and we show that the behavior of the 
high-energy emission of this burst can be explained by a model in which the prompt emission consists of two 
components: one is the MeV component due to the synchrotron-self-Compton radiation of electrons accelerated 
in the internal shock of the jet and the other is the high-energy component due to inverse Compton scattering 
of the photospheric X-ray emission of the expanding cocoon off the same electrons in the jet. Such an external 
inverse Compton effect could be important for other Fermi-LAT GRBs, including short GRBs as well. In this 
model, the delay timescale is directly linked to the physical properties of GRB progenitor. 



1. Introduction 

Gamma-ray bursts (GRBs) were only sparsely ob- 
served in the > 100 MeV range, until the Fermi satel- 
lite was launched on June 11 2008 [H]. Now Fermi 
provides extremely broad energy coverage, 8 keV — 
300 GeV, with high sensitivity for GRBs, and is ac- 
cumulating a wealth of new data which open a com- 
pletely new window on the physics of GRBs. The 
high-energy temporal and spectral data provided by 
Fermi can severely constrain the physical parameters 
of the GRB emission region and the circumburst en- 
vironment, which will lead to a deeper understanding 
of the central engine and the GRB progenitors, and 
will also constrain models of high-energy cosmic ray 
acceleration 0,0]. 

GRB 080916C has the largest isotropic 7-ray energy 
release so far, E Jt i so ~ 8.8 x 10 54 erg (with redshift z ~ 
4.35). Fermi LAT obtained its high-energy emission 
data quite extensively, showing several important new 
properties Q: 

(i) The time-resolved spectra (with resolution ~ 5— 
50 s) are well fitted by a smoothly broken power- 
law function (the so-called Band function) from 
8 keV up to a photon with energy ss 13.2 GeV. 

(ii) The e > 100 MeV emission is not detected to- 
gether with the first e < 1 MeV pulse and the 
onset of the e > 100 MeV emission coincides 
with the rise of the second pulse (~ 5 s after the 
trigger) . 

(iii) Most of the emission in the second pulse shifts 
towards later times as higher energies are con- 
sidered. 



(iv) The e > 100 MeV emission lasts at least 1400 s, 
while photons with e < 100 MeV are not de- 
tected past 200 s. 

Some other Fermi-LAT GRBs also display high- 
energy lags, similar to the properties (ii) and/or (iii) 
[1, [J, and then they should be very important to un- 
derstand the prompt emission mechanism of GRBs. 
We will call the e % 1 MeV emission and the e > 
100 MeV emission "MeV emission" and "high-energy 
emission" , respectively. 

A simple physical picture for the property (i) is that 
the prompt emission consists of a single emission com- 
ponent, such as synchrotron radiation of electrons ac- 
celerated in internal shocks of a relativistic jet. In 
this picture, the peak of the MeV pulse could be at- 
tributed to the cessation of the emission production 
(i.e., the shock crossing of the shell) and the subse- 
quent emission could come from the high latitude re- 
gions of the shell. Thus the observed high-energy lag 
for the second pulse (property (iii)) requires that the 
electron energy spectrum should be harder systemat- 
ically in the higher latitude region. This would imply 
that the particle acceleration process should definitely 
depend on the global parameters of the jet, e.g., the 
angle-dependent relative Lorentz factor of the collid- 
ing shells, but such a theory has not been formulated 
yet. The property (ii) could be just due to the fact 
that the two pulses originate in two internal shocks 
with different physical conditions for which the elec- 
tron energy spectrum of the second internal shock is 
harder than that of the first one [H . 

Another picture is that the prompt emission con- 
sists of the MeV component and a delayed high-energy 
component. The latter component could be produced 
by hadronic effects (i.e., photo-pion process and pro- 
ton synchrotron emission) [|| 0] , but they require ex- 
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tremely large total energy budget Q. 

In this paper, we discuss a different two-component 
emission picture in which the delayed high-energy 
component is produced by leptonic process (i.e., elec- 
tron inverse Compton scattering). We focus on the 
effect that the ambient radiation up-scattered by the 
accelerated electrons in the jet can have a later peak 
than that of the synchrotron and synchrotron-self- 
Compton (SSC) emission of the same electrons (cor- 
responding to the property (hi)) [t| QUEH- Provided 
that the seed photons for the Compton scattering 
come from the region behind the electron accelera- 
tion region of the jet (see Figure [TJ, the up-scattered 
high-energy photon held is highly anisotropic in the 
comoving frame of the jet, i.e., the emissivity is much 
larger for the head-on collisions of the electrons and 
the seed photons. As a result, a stronger emission is 
observed from the higher latitude regions, and thus 
its flux peak lags behind the synchrotron and SSC 
emission. 

Here we propose that the seed photons may be pro- 
vided by the photospheric emission of an expanding 
cocoon. GRB 080916C is a long GRB, and it may 
originate from the collapse of a massive star. The rel- 
ativistic jet produced in the central region penetrates 
the star and deposits most of its energy output into 
a thermal bubble, or cocoon, until it breaks out of 
the star [13, EH EH EEl- The cocoon can store an 
energy comparable or larger than the energy of the 
prompt emission of the jet, and thus it may make an 
observable signature outside the star [TtJ ■ The co- 
coon escaping from the star will emit soft X-rays, and 
these can be up-scattered by the accelerated electrons 
in the jet into the high-energy range. The optical thin- 
ning of the expanding cocoon may be delayed behind 
the prompt emission of the jet, so that the onset of 
the high-energy emission is delayed behind the MeV 
emission (corresponding to the property (ii)). Thus 
this model has the potential for explaining the two 
delay timescales; the delayed onset of the high-energy 
photons (property (ii)) is due to the delayed emission 
of the cocoon, while the high-energy lag within the 
second pulse (property (hi)) is due to the anisotropic 
inverse Compton scattering. We also show that the 
combination of the time-averaged spectra of the SSC 
and the up-scattered cocoon (UC) emission is roughly 
consistent with the observed smooth power-law spec- 
trum (property (i)) (see Figure [5]). 

As we will explain below, the UC emission is short- 
lived and may not account for the whole high-energy 
emission, which lasts much longer than the MeV emis- 
sion (property (iv)). It is natural that the high-energy 
emission in later times is related to the afterglow, i.e., 
produced by the external shock which propagates in 
the circumburst medium. This possibility is stud- 
ied by [HI]. They claim that even the onset phase 
of the high-energy emission is produced by the ex- 
ternal shock. However, the rise of the flux of GRB 
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Figure 1: Geometry of our model. The jet with opening 
angle 6j ~ 0.01 and the cocoon with opening angle 
6 C ~ 0.8 are ejected from the collapsing star. The cocoon 
becomes optically thin at r = r p h some time (<• 5 s) after 
the burst trigger. The second pulse of the prompt 
emission is produced by the accelerated electrons in the 
internal shock of the jet at r — Vi, and the cocoon X-ray 
photons are up-scattered by the same electrons at r — n 
into the GeV range. 

080916C in the LAT energy range (~ t 6 ) is too steep 
for the external shock to reproduce it. Thus at least 
the first part of the delayed high-energy component of 
this burst should be related to the prompt emission. 

Although we focus on GRB 0809 16C to examine 
whether the UC emission is viable for its properties in 
the high energy range in this paper, our modeling is 
generic, and the UC emission could be important for 
some other Fermi-L AT GRBs (including short GRBs). 
We will discuss the temporal and spectral properties 
of the cocoon emission (in § 2), the synchrotron and 
SSC emission of the jet (in § 3), and the UP emission 
of the jet (in § 4). We summarize our model and its 
implications in § 5. For more details, see [l9| . 



2. Cocoon Emission 

We take the total energy and the total mass of 
the cocoon and the stellar radius as parameters, E e , 
M c , and r*, respectively. After the jet breaks out of 
the star, the cocoon expands with the sound speed, 
c s = c/\/3, where c is the speed of light. The cocoon 
expands by its own thermal pressure in the comoving 
frame as expected in the standard fireball model 0]. 
Its expansion speed in the comoving frame suddenly 
becomes close to c, and then the opening angle of the 
cocoon measured in the central engine frame is given 
by 9 C ~ rjg X , where F Ciex is the Lorentz factor cor- 
responding to c s . Thus we obtain C ~ \/2/3 ~ 0.8. 
This rough estimate is consistent with the results of 
the numerical simulations [151 ]. 

The cocoon material accelerates as T cx r and 
reaches the terminal Lorentz factor T c = E c /M c c 2 
at r s = r,r c = 5.0 x 10 12 r»,n(r c /50) cm. (We adopt 
the notation Q = 10 X Q X in cgs units throughout this 
paper.) For r > r s , the cocoon material is dominated 
by the kinetic energy. The radiation stored in the 
cocoon is released when the opacity for the electron 
scattering becomes less than unity. The photosphere 
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radius is given by 
r p h - 



E c a T 



2w(l — cos 6 c )T c m p c 2 
■50 



1/2 



2.1 x 10 14 E l c % 



■1/2 



cm, 



(1) 



where <jt is Thomson cross section. 

The cocoon may become optically thin later than 
the onset of the first MeV pulse (i.e., the burst trigger 
time). Let t = be the photon arrival time at the 
earth emitted at the stellar surface at the jet breakout. 
The first MeV pulse produced within the jet at r = 
Ti is observed at t ~ Atj = rj(l + z)/(2cr 2 ), where 
Tj is the bulk Lorentz factor of the jet and z is the 
source redshift. This timescale is comparable to the 
angular spreading timescale of the pulse, and we can 
take Ati ~ 2 s for GRB 080916C. The second MeV 
pulse is observed ~ 5 s after the burst trigger, i.e., at 
f ~ At* + 5 — 7 s. On the other hand, the cocoon 
photospheric emission is observed at 



t~ At c 



2cYl 



(1 + 2) ~ 7.5 E, 



1/2 
c,52 



-5/2 



(2) 



where we have used z ~ 4.35. The cocoon photo- 
spheric emission may be observed from t = At c to 
t = At c + Atd, where Atd ~ At c is the time during 
which the cocoon will be adiabatically cooled. If in- 
ternal dissipation occurs in the jet at r = > r p h 
making the second MeV pulse within the duration of 
the cocoon emission, the cocoon photons may be up- 
scattered to higher energies by the energetic electrons 
within the dissipation region of the jet, which may be 
observed along with the second MeV pulse. Therefore 
we require a condition Ati < At c < Ati + 5. This con- 
dition puts a constraint to the physical parameters of 
the cocoon, 



0.3 < E. 



1/2 

c,52 



r \ -5/2 



(3) 



We adopt the parameters E c ~ 10 52 erg and r c ~ 50 
for the purposes of calculating the flux of the cocoon 
emission. 

The comoving temperature of the cocoon when 
its opening angle becomes 6 C is approximately given 
by T( Dit ~ [£ c /(27r(l - cosf^a)] 1 / 4 , where a 
is the Stefan constant. Then the temperature 
at the photosphere is estimated to be T^ h = 
T i ' nit (r s /r >t ) _1 (rph/r ;j ) _2 / 3 . Non-thermal electrons in- 
jected by internal shocks within the photosphere may 
make the emission be quasi-thermal. Its spectum is 
written by 




for e < e 



for e;i 



ph 

< e < e\ 



(4) 



CO 

cut ) 



where and F£° h are given by 



2.*> kr; h ^- ~ 1.2 e^t^ 1 k,\" ( - 



(5) 



f: 



{l + zf2 ^Hf f rph 

« K1 ph l ' 



- 31 K% r *,ii keV cm ~ 2 s_1 keV" 1 , (6) 

where — s^/h (h is the Planck constant), and we 
have taken the luminosity distance of GRB 080916C 
as oIl ~ 1.2 x 10 29 cm. Some numerical calculations 
of the radiative processes in the cocoon show f3 ~ — 1 
and e» ~ 30 x e c p ° O . 

The observation of GRB 080916C shows that there 
is no excess from the Band spectrum at the X- 
ray band, > 10 keV, and we obtain a rough up- 
per limit of the cocoon X-ray emission £p\F^° b < 



40 keV cm" 



This limit leads to another con- 



straint on the cocoon parameters, 



r*,n > 0.8 E* ( ^ 



(7) 



3. Synchrotron and SSC Emission 

We can constrain, from the Fermi observation, the 
global physical parameters of the jet: the bulk Lorentz 
factor Tj , the emission radius of the second pulse r, , 
and the opening angle 9j. First of all, from the ab- 
sence of a 77 absorption cutoff, we obtain a lower 
limit on the bulk Lorentz factor of the jet, Tj > 870 
Q • Since the angular spreading timescale of the pulse 
is Ati — 2s for the second pulse, similar to the first 
pulse, the emission radius is estimated by 



2.2 x io 16 r 2 3 



At\ 

- cm. (8) 



Since GRB 0809 16C is so bright, it is probable that 
the jet is viewed on-axis, and we adopt this assump- 
tion as a simplification (see Figure[T]). In this case, the 
cocoon is viewed off-axis, since the jet cone will not 
be filled with the cocoon material. The cocoon emis- 
sion is thus less beamed, but this off-axis effect is not 
significant because the opening angle of the jet can be 
estimated to be small. The isotropic 7-ray energy of 
this burst is 8.8 x 10 54 erg. To obtain a realistic value 
of the collimation-corrected 7-ray energy, < 10 51 erg, 
the jet opening angle is constrained by 6j ^ 0.015. We 
adopt 0j ~ 0.01, and having adopted a nominal value 
of T c ~ 50 in accord with the observed time delay of 
the high-energy emission (equation [3]), we obtain 



0.5 < 1. 



(9) 
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Thus the off-axis dimming and softening effects are 
not significant for the cocoon emission. 

We assume that the jet is dominated by the ki- 
netic energy of protons and we estimate the physi- 
cal parameters of the jet dissipation region for the 
second MeV pulse. In the scenario where the dissi- 
pation is due to internal shocks [2], the collisionless 
shock waves can amplify the magnetic field and ac- 
celerate electrons to a power-law energy distribution, 
which then produce synchrotron radiation and SSC 
radiation. At r — r, ; , the comoving number density of 
the jet is estimated by n' = L/Stij (Aitr\ m p c 2 (l + z)), 
where L is the isotropic-equivalent luminosity of the 
jet. The internal energy density produced by the in- 
ternal shock is given by u' = n'O p m p c 2 , where 9 P is 
a factor of order unity. Assuming that a fraction es 
of the internal energy of the protons is carried into 
the magnetic field, the field strength is estimated by 
B' = (Sttcbu') 1 / 2 . Assuming that a fraction e e of the 
internal energy of the protons is given to the electrons, 
the minimum Lorentz factor of the electrons is given 
by 7m = [ip - 2)/(p - l)](m p /m e )e e O p , where p is 
the index of the electron energy distribution and we 
have assumed that it is similar to the value inferred 
from the spectrum of the first pulse, p ~ 3.2. We 
find that the cooling electron Lorentz factor 7 C ~ 7 m 
for fitting the observational data. In this case we ob- 
tain 7c ~ 360(e B /10^ 5 )- 1 / 3 (r/4 x 1CT 4 )- 2 / 3 , where 
t = arn'ri/Tj is the Thomson optical depth. 

The synchrotron characteristic energy and the syn- 
chrotron peak flux (at the synchrotron energy corre- 
sponding to 7 C ) are estimated by 



3heB' 2 2Tj 

47TTO e C 7m 1 + 2 



2 7 L 1/2 r- 2 



,6/a ( JflV'VieV eV 
p Vio-v \0AJ ' 



(10) 



y^gr^) ^ i3xio4 , 

2 47r ^2 55 j,3 



m e c 



(yjpg) keV cm' 2 s' 1 keV" 1 (11) 

where N = [LAt l /(l + z)]/(I\,to p c 2 ). The lst-order 
SSC characteristic energy and the SSC peak flux are 
approximately 



l/2 r -2 



-fn — ^lm £ m — 1-7 £55 1^3 { } ^ 



s ¥ 2 & 2 (^y^, (i 2 , 



We find that the lst-order SSC radiation can account 
for the observed MeV emission of this burst. The 
emission just below e c is suppressed by synchrotron 
self-absorption effect, and the 2nd-order SSC emission 
is suppressed by the Klein-Nishina effect. 



4. Upscattered Cocoon (UC) emission 

Here we derive the observed spectrum of the UC 
emission as a function of the polar angle 9 of the emit- 
ting region on the shell (see Figure [T]). The spectrum 
of radiation scattered at an angle 0' sc relative to the di- 
rection of the photon beam in the Thomson scattering 
regime is given by [|J [ljj, HTj] 

j'AKc) = ^t(I-cosC) / d 7 N'( 7 ) f dyJ £ ,(l-2y+2y 2 ), 



(14) 

where y = e'/[2j 2 e' s (l — cos6' sc )]. This is the scat- 
tered radiation emissivity in the jet comoving frame, 
AT' (7) is the electron energy spectrum, and J' e , is the 
intensity of the seed photons averaged over solid an- 
gle, i.e., the mean intensity. The term £ = 1 — cos6' sc 
describes the anisotropy of the spectrum, and this is 
due to the fact that the IC scattering is strongest for 
the head-on collisions between electrons and seed pho- 
tons. This implies that the UC emission in the ob- 
server frame is stronger from the high-latitude region 
of the shell, so that its flux peak lags behind the on- 
set of the synchrotron and SSC emission of the same 
electrons, which have isotropic energy distribution in 
the comoving frame of the jet. 

In order to concentrate on the time-averaged spec- 
trum including the high-latitude emission, we calcu- 
late the flux of the UC emission by neglecting the 
radial structure of the emitting shell for simplicity. 
The peak energy and the peak flux of the UC emis- 
sion are reduced as functions of the angle parameter 
q{6) = T 2 6 2 : 



-UC 



27^ m 



1 + T 2 9 2 



160 MeV 



4q 



Q_„, „, _CO 77CO 

oT7 m 7 c e ph ^ £ph ^ ^ p202p 



580 keV cm" 2 s" 1 



40g 2 



V400/ V400/ V 30 keVcm^s- 1 



(1 + <Z) £ 

t ph r e ph 



4 x 10- 4 / 

(16) 



f? c ~ tF £ ~ 3.4 L^ 2 rr 8 



J 55 1 j,3 I 2 
x el' 2 (y^s) 172 keV cm- 2 s" 1 keV' 1 ^) 



where the functions in the brackets [ ] both have val- 
ues of zero at q — and q — 00 and have peaks of 1 at 
q = 1 and ~ 1.4 at q = 2/3, respectively. This means 
that the UC flux has a peak at q m 1, or 6 ~ 
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Figure 2: Time-averaged spectrum of the second pulse 
calculated in the up-scattered cocoon emission model. 
The lst-order SSC component plus 2nd-order SSC 
component without including the Klein-Nishina effect 
(dashed line), the cocoon photospheric emission 
(dot-dashed line), and the UC emission (dotted line) are 
shown. The thick solid line represents the combination of 
these, taking account of the Klein-Nishina effect, which is 
roughly consistent, at e > 1 MeV, with the Band model 
spectrum (thin solid line) with 95% confidence errors 
(dot-short-dashed lines) (from the LAT/GBM group of 
Fermi). The bump at ~ 30 GeV is so dim as not to be 
detected. The adopted parameters are listed in equation 
QUI and (fT8l). 

i.e., the peak time of the UC emission lags behind 
that of the SSC emission on the angular spreading 
timescale, At; ~ 2 s. This is consistent with the ob- 
served lag of the GeV emission onset behind the MeV 
emission peak of the second pulse of GRB 080916C. 
Here the values of the jet parameters t = 4 x 10~ 4 and 
7m = 7c = 400 are applicable for the lst-order SSC 
emission of the jet being consistent with the observed 
MeV emission component (see § 2). This indicates 
that the UC emission of the electrons accelerated in 
the internal shock of the jet, emitting the observed 
MeV emission, can naturally explain the observed flux 
in the GeV range. 

If the flux of the cocoon photospheric X-rays is 
given, i.e., E c , T c , and r*, are given, the fluxes of 
the UC and SSC emission of the jet are determined 
by the jet parameters L,Tj, AU,eB, and e e . Since 
Ati ~ 2 s is roughly given by the observations, and 
this value is necessary to explain the observed high- 
energy lag timescale, we have four free parameters. 
On the other hand, we have four characteristic ob- 
servables; the peak fluxes and peak photon energies of 
the SSC component and the UC component. There- 
fore the jet parameters are expected to be constrained 
tightly. 

Figure [2] shows the result of the time-averaged spec- 
trum of the second pulse for the cocoon parameters 

£ c ,52 = 1.0, T c = 52, r*,u=2.5, (17) 

and (3 = —1.2. These values satisfy the constraints on 
the cocoon parameters, equations © and ([7]). This 
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Figure 3: Photon number fluxes of the second pulse in 
several frequency ranges calculated in the up-scattered 
cocoon emission model. Each lightcurve is normalized to 
a peak flux of unity. The peak of the GeV lightcurve is 
delayed behind that of the MeV lightcurve. 

figure shows that our model is roughly consistent with 
the observed spectrum at e > 1 MeV. The adopted 
values of the jet parameters are 

L 55 = l.l, I\j )3 = 0.93, Ai, = 2.3s, e B = 10~ 5 , e e = 

(18) 

and p = 3.2. The corresponding values of the opti- 
cal depth for electron scattering and the characteris- 
tic electron Lorentz factors are r = 3.5 x 10~ 4 , j m = 
400, 7 C = 390. Figure [3] shows the results of the 
multi-band lightcurves for the same parameters. Each 
lightcurve is normalized to a peak flux of unity. This 
clearly displays the lag of the high-energy emission 
peak. 

5. Discussion 

We have discussed a model in which the prompt 
emission spectrum consists of an MeV component pro- 
duced by the SSC emission of electrons accelerated in 
internal shocks in the jet and the high-energy compo- 
nent produced by up-scattering of the cocoon X-rays 
off the same electrons (UC emission), and we have 
shown that this model can explain the above three ob- 
served properties (i), (ii), and (iii) listed in § 1. The 
expanding cocoon may become optically thin some 
time later than the first internal shock of the jet (equa- 
tion [2]), so that the first MeV pulse may not be asso- 
ciated with the UC emission while the second MeV 
pulse may be associated with it (property (ii)). The 
UC emission has an anisotropic energy distribution in 
the comoving frame of the jet so that the observed UC 
emission is stronger from the higher-latitude region of 
the shell. This results in the lag of the flux peak of the 
UC emission behind the MeV emission onset on the 
angular spreading timescale (property(iii)). Figure [2] 
shows that the combination of the SSC and UC emis- 
sion can reproduce the observed high-energy spectral 
data (property (i)). The UC emission is short-lived 
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(roughly for At c ~ 7.5 s) and may not account for the 
whole high-energy emission which lasts longer than 
the MeV emission (property (iv)). It is natural that 
the high-energy emission in the later times is related to 
the afterglow. This has been shown by [l8j]. However, 
the early portion of the high-energy emission should 
be the UC emission, because the external shock can- 
not reproduce the observed steep rise of the flux. 

We have focused on GRB 080916C because of its 
extensive data on the high-energy emission. How- 
ever our model is generic, and it could apply to 
other Fermi-LAT GRBs with typical parameters L ~ 
10 53 erg s"\ 6j ~ 0.1, and Tj ~ 300 for which the 
prompt emission in the soft 7-ray range is produced 
by the lst-order SSC radiation of electrons [13]. Some 
short GRBs might originate from the collapses of the 
massive stars [21|, and thus they could have the de- 
layed UC emission. Even if other short GRBs are pro- 
duced by the compact star mergers, it might be pos- 
sible that the jet is accompanied by the delayed disc 
wind [22j and that the emission from the disc wind 
is up-scattered by the electrons accelerated in the jet. 
For either progenitor models, the delayed high-energy 
emission associated with short GRBs, if any, would 
provide an interesting tool to approach their origins. 

A simple prediction of our model is that prompt 
emission spectra of some GRBs would have an excess 
above the Band spectrum around ~ 1 keV due to the 
cocoon photospheric emission, and this excess should 
have a different temporal behavior from that of the 
MeV emission. In addition, we expect GRB 080916C 
to have had bright synchrotron emission in the optical 
band, like the "naked-eye" GRB 080319B. 

If our model is correct, we can constrain the pa- 
rameter range for which hadronic effects are impor- 
tant on the high-energy emission of GRBs, and we 
can also constrain the models of high-energy cosmic 
ray acceleration. Also, the delay time of the onset 
of the high-energy emission is directly linked to the 
optical-thinning time of the expanding cocoon, which 
constrains the physical parameters of the progenitor 
star and the cocoon material of GRBs. For GRB 
080916C, the stellar radius r* and the total energy 
E r and mass M r of the cocoon are constrained to 



> 



be 0.3 < E-' 2 (M c /10- 4 M Q y- 5 < 0.9 and r* >n 
0.8 El 52 (M c /10- 4 M Q )- 3 (see equations 1 and EJ. 
The cocoon energy and the cocoon mass come from 
the jet energy released within the star and the stellar 
mass swept by the jet, respectively. These constraints 
therefore provide potential tools for investigating the 
structure of the progenitor star just before the ex- 
plosion, as well as the physical conditions of the jet 
propagating inside the stellar envelope through either 
analytical [12, EH or numerical 0, EB] approaches. 
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